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THE  FORCED  OSCILLATION  OF  SHALLOW  DRAFT  SHIPS 
By  W.  D.  Klffl' 

Considering  a  threa-diasnsional  disk  with  a  circular  or 
elliptic  waterplane  as  a  ship  of  shallow  draft,  the  interaction  between 
the  disk  oscillating  with  six  degrees  of  freedom  and  the  induced  motion 
of  the  surrounding  fltiid  is  investigated*  The  potential  problem  is 
formulated  by  the  use  of  the  concept  of  surface  distributed  sources 
so  that  an  integral  equation  for  the  source  density  is  obtained*  The 
nmerleal  solution  of  the  integral  equation  is  found  and  the  dependence 
of  added  mass,  added  moment  of  inertia,  and  damping  factors  on  various 
frequencies  is  determined* 

I*  INTRODUCTION 

Consider  the  forced  oscillation  of  a  rigid  disk  on  a  free 
surface*  The  form  of  the  disk  is  such  that  it  will  create  large  wave* 
making  effects*  The  present  analysis  thus  serves  as  a  complement  to 
the  thin-ship  theory  Cl,  33^*  When  the  draft  is  small,  the  normal 
derivative  of  Qreen's  function  (the  potential  due  to  a  unit  source) 
does  not  enter  into  the  kernel  of  the  Integral  equation,  and  moreover 
the  kernel  can  be  evaluated  e^llcltly*  However,  for  arbitrary  cross- 
sections  C2],  the  presence  of  the  normal  derivative  of  Green's 


^  Staff  Member  of  the  Flight  Sciences  Laboratory,  Boeing  Scientific 
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*  Numbers  in  brackets  designate  References  at  end  of  paper* 


function  mnkca  the  numerical  solution  extremely  difficult* 

It  is  assumed  that  the  amplitude  of  oscillation  is  sufficiently 
small  so  that  the  problem  can  be  linearized*  One  may  choose  the  ratio 
of  actual  draft  to  the  wavelength  c  as  a  small-parameter  and  develop 
a  perturbation  procedure  following  the  method  of  MacCamy  [33*  It 
then  turns  out  that  the  added  mMS,  added  moment  of  lnertia«  and  damp¬ 
ing  factors  associated  with  surge*  sway*  and  yaw  are  terms  of  order  c* 
Furthermore*  this  procedure  reduces  the  potential  problem  of  an  oscillat¬ 
ing  disk  to  a  series  of  boundary  value  problems  of  increasing  orders 
in  c* 

The  solution  of  the  lowest  order  problem  is  sought  in  the 
form  of  a  potential  due  to  the  distribution  of  sources  over  the  disk* 
Then*  the  singular  property  of  surface  distributed  sources  enables  one 
to  determine  the  density  from  an  integral  equation  which  satisfies  the 
boundary  values  on  the  disk*  The  numerical  solution  of  this  integral 
equation  can  be  obtained  approximately  by  use  of  a  finite  set  of  linear 
equations  which  relate  the  unknown  densities  at  selected  pivotal  points* 

A  scheme  has  been  developed  for  evaluating  the  singular  inte¬ 
gral  in  two  variables*  (The  fundamental  singularities  are  log  r  and  •*) 
From  the  solutions  of  the  linear  equations*  pivotal  values  of  pressure 
can  be  obtained*  Thus*  knowing  the  pressure  distribution*  one  can 
evaluate  added  mass*  added  moment  of  inertia*  and  damping  factors  for 
heave*  roll,  and  pitch  corresponding  to  a  specific  frequency*  It  was 
found  that  accurate  results  for  high  frequency  notion  require  a  large 
number  of  pivotal  points* 


II.  QgNERAL  FORMnLATIOK 


Suppose  the  plane  y  «  0  represents  the  undisturbed  free  sur¬ 
face  of  an  inTlseld  amd  incompressible  fluid;  0  x  and  0  z  are  axes  which 
lie  in  that  surface*  The  hull  of  a  ship  in  the  equilibrium  position 
is  given  by 

y  *  S®(x,  z)  ,  (2.1) 

where  the  form  of  the  ship  is  assumed  to  satisfy  the  following  condi¬ 
tions: 

S^(Xf  z)  B  0  along  the  edge  of  the  ship,  j 

S^(x,  z)  ■  c,  maxlmtim  draft  at  the  center,  >  (2.2) 

S°(x,  z)  ■  S®(-x,  z),  and  S°(x,  z)  *  S®(x,  -z).  ) 

The  last  two  conditions  Indicate  that  in  the  equilibrium  position,  the 
center  of  the  waterplane  coincides  with  the  origin  0  of  a  time-indepen¬ 
dent  reference  frame  which  is  called  the  space  frame  0  i^,z  * 

Now,  suppose  the  ship  is  excited  into  small  oscillations* 

After  transient  disturbances  have  been  dissipated,  the  resulting  fluid 
motion  will  be  periodic  with  a  frequency. of  oscillation  o*  Considering 
such  a  motion  to  be  irrotatlonal,  the  velocity  potential  may  be  written 
in  the  fora 

0(x,y,z;t)  ■  Re[U(x,y,z)e“^®*]  .  (2.3) 

At  any  Instant,  the  position  of  the  oscillating  ship  may  be  described 
in  terms  of  a  position  vector  JR  ■  iX  jY  'f  kZ,  and  the  Eulerian  angles 
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Ob  10  ^jO  -fkO,  wh«r«  1,  jt  ^  d«not«  unit  vectors  along  axes  of 
“•  X  y  B 

the  space  frame*  One  calls  the  components  X,  f,  and  Z,  surge,  heave 

and  sway,  respectively,  and  the  angular  components  0  ,  0  ,  and  0  , 

X  y  z 

roll,  yaw  and  pitch,  respectively*  Since  the  motion  of  the  ship  Is 
periodic  In  time,  one  has 

R(t)  «  Re[R°e“^®^] 

©(t)  »  RefeV^®^] 

where  and  9^are  euaplltudes  of  the  linear  and  angular  displacements, 
respectively*  From  rigid  body  dynamics,  the  velocity  of  a  particle  on 
a  ship  is  given  by 

£.  ~  R  *  ^  ^  L,'  ♦  (2*5) 

where  jr  Is  the  position  vector  with  respect  to  the  space  frame,  while 
Is  the  position  vector  with  respect  to  a  moving  frame  attached  to 
the  ship* 

In  order  to  show  the  physical  significance  of  small  oscllla* 
tloxis.  Independent  and  dependent  variables  are  made  dimensionless  by 
dividing  each  variable  by  the  wave  length*  The  space  coordinates  then 
become 


;1 


(2*4) 


XBkx,  y“ky»  8  =  k*,  (2*6) 

where  k  represents  the  wave  number  which  Is  equal  to  o^/g  (or  2tT/X), 
g  being  the  acceleration  of  gravity,  and  X  being  the  wave  length  of 
free  waves  of  frequency  o*  Assuming  the  ratio  of  actual  amplitudes  to 
the  wavelength  Is  small,  one  may  express  the  normalized  amplitudes  for 
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th«  linear  displaeements  using  a  small  parameter  a  as  ■  aX^, 

T**  ■  aT^f  and  ■  aZ^.  Similarly,  the  amplitudes  for  small  angular 
displacements  may  be  expressed  as  0^  m  0^  *  a0^|  and  0^  ■  a0^# 

X  X  y  jT  Z  Z 

2 

Accordingly,  neglecting  the  terms  involring  a  ,  the  velocity  potential 
for  the  resulting  fluid  motion  can  be  expressed  as 

0  ■  a0^  ,  (2.7) 

where  0^  and  its  derivative  are  bounded.  I'urthermore ,  the  linearised 
dynamic  pressure  is 

TT  ■  -  p0^(x,y,s;t)  ■  -  »  (2.8) 

where  p  is  the  density  of  the  fluid,  and  the  linearised  surface  eleva¬ 
tion  becomes 

n  ~  i  0^(;,O,S;t)  »  -  I  0j(J.O,f;t)  .  (2.9) 

Since  the  fluid  is  considered  to  be  Incompreseible,  the  velocity  poten¬ 
tial  must  satisfy  the  Laplace  equation, 

7^0  «  0  , 

or 

.  0  in  y  <  0,  outside  S  .  (2.10) 

At  a  free  surface,  the  first-order  boundary  condition  to  be  satisfied 

is 

0^t  ♦  •  0  . 

or 

0-  - 
7 


kU  «  0 


on  y  s  0,  outside  S 


(2.11) 
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Next  I  on  the  Immersed  surface  the  kinematic  condition  can  be  written 


0-  =  £(x,y,a!t)  Hx,y,£:t)  on  S(x,z:t),  (2*12) 


where  jn  represents  the  unit  vector  normal  to  the  surface,  l*e*, 
n  e  in^  +  jOy  ♦  velocity  vector  of  a  point  on  the 

surface*  Note  that  the  surface  S  of  the  ship  in  motion  differs  from 
the  surface  at  rest  which  was  Introduced  earlier*  Since 


[0-^  *  ft  ^  » 


it  is  consistent  with  the  omission  of  the  second  order  terms  to  require 
that  condition  (2*12)  be  satisfied  on  S°*  By  substituting  (2*5)  into 
the  right  hand  side  of  (2*12)  one  obtains 

e  s  e 

0*KXa  'fTn  ^  Zn  •fdq  -fOq  4-6q 
X  y  z  x^x  y^y  z^z 


E  [U  (x.y,®)]-  = 
j«l  ^  “ 


-io[X®n  ♦  Y**n  ♦  Z®n  +  0®q  ♦  <l_  1 

X  y  z  X  X  y  j  z  z 


on  S®(x,5)  * 


(2*13) 


Since  the  problem  is  linearized,  the  total  potential  U  is  ejqpressed  as 

6 

E  U .  J  s  1,  2,  —  6,  and  (q  ,  q„,  q.)  represents  the  components  of 
j=l  J  X  y  z 

a  vector  given  by  ^  £  x  ^* 

Finally,  at  large  distances,  the  propagating  disturbance  must  have  the 
form  of  a  radically  outgoing  wave,  l*e*. 


.  XiSi  .''J  .  (A) 


(2*14) 
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wher*  r  ■  Jx^  *  and  0  «  arc  tan  . 


To  show  clearly  the  dependence  of  the  solution  on  paraaeters 
describing  the  type  of  oscillation,  one  relates  the  dynamic  pressure 
functions  Pj(x,y,s)  j  ■  1,  2,  — *  6,  to  the  potential  function 
n(x,y,i)  by 

rx®  Y®  Z®  ^  cO  *■  efi  *  ft®  1  <  ^  r  Ti  Z 

«LT  *—  'z  *—  »j  V4  *k  Vs  VsJ  ■  • 

(2.15) 

Than  the  boundary  conditions  (2.10),  (2.11),  (2.13),  and  (2.14)  can  be 
transformed  Into 


:  0 

in 

y 

<  0  ,  outside  S 

(2.16) 

'  • 

0 

on 

y 

B  0  ,  outside  S 

(2.17) 

''I’n 

(Pz'n  * 

“y  ’ 

‘^S’n  ■  \ 

<P4>n 

z 

-I“y  * 

Cp,), 

3 

1 

Z  X  [ 

1  “x  "  a  “»  »  > 

(2.18) 

-in 

a  X 

on 

S®(x,*)  ,  / 

Pj  -  «  r  -  -  . 


(2.19) 


where 


/1^(S^)S(SJ‘ 


'  X  X 


,  and  n  ■ 


L+(S  )^+(S,)^ 

X  % 


(2.20) 


III.  SHALLOW  DRAFT  APPROXIMATION 


In  this  section,  the  first  order  forces  and  moments  due  to 
the  dynamic  pressure  are  related  to  the  added  mass,  added  moment  of 
inertia,  and  damping  factors  of  an  oscillating  ship.  The  perturba¬ 
tion  procedure  for  a  ship  of  shallow  draft  is  developed. 

The  forces  and  moments  due  to  the  dynamic  pressure  exerted 
on  the  immersed  surface  are  given  by 


F  =  JJ  TT(x,y,z;t)»n  dS  ,  (J.l) 

S(x,z;t) 

and 

G  s  JJ  TT(x,y,z;t)»^  dS  .  (3*2) 

S(x,z; t) 


Since,  from  (2.8),  TT(x,y,z;t)  «  -  ap0^,  auid  since  S(x,z:t)  ■  S^-foS^, 
the  first  order  forces  and  moments  are 


£■  n  -  ap  JJ  0j(;,y,;;t)«n  dS  ,  (3.3) 

S®(x,z) 

and 

0^  «  -  ap  JJ  0j(x,y,z;t).^  dS  .  (3«^) 

S®(x,z) 


9y  (2.3)  mnd  (2.13),  one  finds 

0j  «  -  g  Re[R*-pje"^®*3,  where  j  »  1,2,3,  (3*5) 

or 

e  -  g « ReCS^Pje"^*^^  ],  where  j  ■  4,3,6. 


(3.6) 
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Therefor*,  the  substitutions  of  (3*9)  into  (3*3),  nnd  (3*6)  into 
(3*4)  yield 

£*•  »  apg  JJ  R*[R^Pje“^®^]*n  dS  ,  where  J  -  1,2,3,  (3»7) 

S®(x,z) 

and 

s  apga  JJ  ReC^p  .e”^®^]*jq  dS  ,  where  J  ■  4,3,6*  (3*8) 

S®(x,a) 

Suppose  the  first  order  forces  and  mofflents  are  ejqpressed  in 

the  fora 

.  <3-9) 

and 

a^«-YaS-Ha9  ,  (3*10) 

where  M  and  ^  are  called  the  added  mass,  and  added  moment  of  Inertia, 
while  ^  and  ^  are  called  the  damping  coefficients  for  translation  and 
rotation,  respectively* 

^y  substituting  (2*4)  into  (3*9)  and  (3*10)  one  obtains 

»  o  R*C(o^  +  ioN)^*-*"^®^  ]  ,  (3*11) 

and 

0^  »  a  R*C(o^  ♦  ioH);  .  (3*12) 

In  equating  (3*7)  to  (3*11),  and  (3*8)  to  (3*12),  all  variables  are 
made  dimensionless  by  the  transformation  x  *  kx,  y  ■  ky,  and  s  «  ki. 


then  one  has 
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■  JJ  Pj*ii  k^I^po  ■  I®{JJ  Pj*E  *^®}*  *li®**«  J  ■  1«2,5, 

S®(x,z)  S®(x,z) 


(3.13) 


k'*I/p  .  8.{  JJ 


Pj’iL 

(x,z) 


^“{JJ  Pi*a'is},  where  j  »  ^.5»6« 
S°(x,z) 


(3.14) 


S^(x,z)  represents  the  image  of  S^(x,z),  and  hence  S^(x,z)  >  k^S^(xtz). 

Next,  it  will  be  shown  how  the  shallow  draft  approximation 
cam  be  developed.  Here  the  shallow  draft  means  that  the  ratio  of 
actual  draft  e  to  the  wavelength  X  is  small.  Therefore,  one  may  begin 
by  writing  the  hull  of  the  ship  in  the  form 

y  =  «S^(x,z)  -a^x^a  -b^z^b,  (3.15) 


with  e  s  2TTe/X  as  the  small  dimensionless  parameter  appropriate  for  the 
came  of  shallow  draft.  S^(x,z)  satisfies  the  following  conditions 

S^(0,0)  =  1,  S^(x,z)  =  S^{-x,z),  and  S^(x,z)  -  S^(x,-z). (3.16) 

From  (3.15)*  Cl  +  (S  )^  ♦  (S  )^]"^  can  be  expamded  into 

X  7* 

2 

[1  ♦  (S  )2  ♦  (8)2]-^  =  1  ♦  (si)2] 

X  2  c  X  2 


+  ^  (-l)“e2®a^C(S^)2  +  (S^)2]“  , 

ms2 


thus  one  obtains 
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-  [-  «sj^  +  0(c^)]  <^+ Cl  ♦0(«^)] C-  cS^  ♦0(«^)]  ~  ,  (3.17) 

wher*  n  ■  -  cS^  +  0(€^),  a  ■  1  ♦  O(e^),  and  n  «  -  «S^  ♦  O(t^)  (3*18) 
XX  y  z  z 

Daalgnatlng  the  right  hand  member  of  (2«l8)  aa  Hj(x,s) 

J  «  1,2,  6,  theae  quantltlea  are  expanded  Into 

«D 

Rj(x,z)  *  ^  e™H®(x,z)  .  (3*19) 


From  (3*18)  one  now  flnda 

R^(x,z)  *  0,  R2(x,z)  b  1,  R^(x,z)  s 

R^(x,z)  ■  0,  Rg(x,z)  ■  i,  I^(x,z)  ■ 

R^(x,z)  B  R^(x,z) 

R^(x,  z)bO,  ••••  • 


0,  R^(x,z)  •  ~ 
-s];,  R^(x,z)  »  0, 


(3.20) 


Sladlarly,  the  preaaure  would  alao  hare  to  be  developed  in  the  fora 

Pj(x,y,z)  B  ^  «“p“(x,y,z)  j  b  i,2,  -  6  .  (3.21) 

BbO 


Here  one  makea  uae  of  Taylor 'a  theorem  and  wrltea  p®  for  y  «  0 


aa 


I 
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m 

p®[x,  €S^(x,z),a]  »  Y,  p"(x,0,*).  (3*22) 

n*0 


Substitution  of  sxpanslons  (3*17),  (3*21)  and  (3*22)  Into 
the  left  hand  member  of  (2*18)  yield 


.n+1 


[-eS^+0(c2)]  Y  1  'k  P“(x,0,*) 


msO  nsO 


dxdy' 


n  *'J' 


on  m 


+  Cl+0(€^)]^  Y  ^  c^^^Cs^Cx.g)  ]°  Pj(x,0,z)  (3»23) 

«_r>  ,.n  *y 


m«0  nsO 


0»  0» 


,n-fl 


♦  C-eS^^O(  e^)  ]  ^  ^  ^  e®^”[S^(x,*)  -^7— p®(x,0,a). 


fflaO  nsO 


dy*^da  ^ 


Upon  equating  coefficients  of  like  powers  of  c  in  (3*19)  and  (3*23)t 
one  obtains 

[p®(x,0,*)3y  *  R®(x,z)  ,  (3*24) 

Cpj(x,0,s)]y  »  R^(x,z)  + S^[p°(x,0,z) ]^  ♦  S^[p®(x,0,*)]^  -  S^[p®(x,0,*)]yy  , 

and  so  forth*  Therefore,  p®  can  be  determined  recurairely  by  solving 


boundary  value  problems 

find  a  function  w(x,y,z) 

Of  the  following  form* 

such  that 

H 

O 

in  y  <0, 

outside  S 

(3.25) 

O 

U 

1 

s’" 

on  y  s  0, 

outside  S 

(3.26) 

Wy  *  L(x,z) 

on  y  B  0, 

within  S 

(3.27) 

,  .  i<®l  .  0(i) 

/7  *• 

as  r  • 

(3.28) 
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Not*  that  hCxfZ)  is  a  function  whoa*  ▼alu*a  ar*  glwan  by  (3*20)* 

In  avaluatlng  th*  add*d  mass,  add*d  monsnt  of  Insrtla,  and 
damping  eosfficlsnts,  th*  half-lsngth  of  th*  ship  a  app*ars  as  an 
Intagral  limit*  Th*r*for«,  all  Isngths  will  b«  mad*  dlm*nslonl*ss 
by  dividing  th*m  by  a  Instaad  of  th*  wav*l*ngth*  This  chang*  can  b* 
readily  carried  out  by  introducing  a  new  set  of  independent  and  de¬ 
pendent  variables f 

X  »  x/a  ,  y  a  y/a  ,  z  =  z/a  , 
and  (3*29) 

w(x,y,z)  ■  aV(x/a,  y/a,  z/a)  where  a  «  ka  • 

The  boundary  value  problems  can  be  stated  as  follows: 

Find  a  fimctlon  V(x,y,z)  such  that 


o 

n 

> 

in  y 

<  0  , 

outside  S  , 

(3*30) 

o 

N 

> 

1 

on  y 

=  0  , 

outside  S  , 

(3*31) 

Vy  »  h(x,z^ 

w 

on  y 

»  0  , 

within  S  , 

(3*32) 

V  -  Xi®I  e»y*^***  . 

0(i) 

as  r 

(3*33) 

Now  the  order  of  th*  added  mass,  added  moment  of  inertia, 
and  damping  factors  involving  a  specific  type  of  oscillation  will  be 
examined  by  th*  us*  of  expansions  which  resulted  from  th*  shallow 
draft  approximation*  Substitution  of  (3*18),  (3*21)  and  (3*22)  into 
(3*13)  snd  (3*14)  lead  to  the  following  expressions: 
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-  0(«), 

k^  /p  B  R«  {JJ  P2(x,0,a)ds}, 

S° 

k^^/p  »  0(«), 

k**iyp  -  -R*{JJ  **pj(x,0,a)ds}, 

S° 

k^Iy/p  «  0(f), 

k^I^/p  ■  R«{JJ  x»pg(x,0,a)dsj’, 

S® 


k^ ^po  ■  0(c), 

k^y/po  =  I“{JJ  p|(x,0,a)ds]-, 

s 

k^^/po  •  0(f), 

k^H^pa  B  »*P^(x,0,»)dsj^, 

S® 

k'^Hy/po  •  0(f), 

k^H^/pO  a  ^“{JJ  3t*Pg(XtO«*)ds}. 


Hfncf,  in  terms  of  the  new  vaxlables,  one  obtains  the  added  mass  and 
damping  coefficient  for  heave  as 


My  B  My/pa^  a  R*{JJ  V2(x,0,?)ds]-  , 

S 

Ny  B  Ny/pa^o  B  Ib{J^  V2(x,0,a)ds]-  , 

S 


(5.34) 


the  added  moment  of  Inertia  and  damping  coefficient  for  roll  as 

B  a  -  ®*{JJ  *^*V^(x,0,Ddsj>  , 

s  (5.35) 

a  H^pa**o  a  -  , 

s 


and  the  added  moment  of  Inertia  and  damping  coefficient  for  pitch  as 

■  ^•{JJ  x»Vg(x,0,z)ds}  , 

"S  (3.36) 

s  H^/pa^o  B  I“{JJ  x»Vg(x,0,Ddsj’  , 

S 
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wh«r«  ?  r«pr«8ant8  the  image  of  S^(x,z)  under  the  trana formation 
X  ■  x/a,  and  z  =  z/a,  and  which  la  equal  to  S®(x,z).  Furthermore* 
alnoe  Pj(x»y»*)  *  aVj(3c*y,z^*  J  «  2*4*6*  from  (5»20)  and  (3»22)  one 
flnda 

[V2(x*0*z')  ]~  a  1,  [V^(x*0*z')  Ip'  a  Z*  and  CVg(x*0*z^V  =  x  on  S 


(3*37) 


Alao*  from  (3*21)*  one  flnda 


[V.(x*0*z) ]~  -  aV.(x*0*z^  *  0  outalde  S 
J  7  J 


(3.38) 


IV.  INTEGRAL  REPRESENTATION 


An  integral  repreaentatlon  of  the  aolutlon  of  the  boundary 
value  problem  which  la  formulated  in  the  preceding  aectlon  will  be 
dlacuaaed  here.  Conalder  a  region  (a  fluid  apace)  bounded  by  the 
htill  of  a  ahip  and  a  free  aurface.  An  oaclllatlng  ahlp  generatea  wavei 
in  thla  region  which  In  turn  react  ao  aa  to  Influence  the  motion  of 
the  ahip.  A  dlaturbance  created  at  a  point  on  the  hull  la  tranamltted 
to  a  point  (x*y*z)  aome  diatance  away  through  the  fluid.  Aa  a  reault 
it  will  be  ahown  that  for  aome  function  f(x*z)  defined  over  the  hull* 
the  reaulting  potential  in  the  region  may  be  expreaaed  in  the  form* 

V(x*y*z)  *  JJ  f(i;,C)G(x*y*z*5*0*C)d5dC  .  (4.1) 

S 


(Hereinafter*  delete  the  tilde  on  varlablea  for  the  aake  of  convenience.) 
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Fhysloallyt  •xpresslon  (4*1)  may  b«  Interpratad  as  tbs  field  poten¬ 
tial  caused  by  distributed  sources  of  the  density  f(x,z)  over  the 
hull  S,  and  then  Green's  function  Q  represents  the  field  potential 
caused  by  a  unit  point  source  at  ((tO,C)  on  the  hull*  The  appro¬ 
priate  Green's  function  [l]  for  the  present  problem  is  obtained  by 
setting  T|  3  0: 

«D 

0(x,y,z»^,0,C)  =  “I  +  ♦ 

0 

•  (4.2) 

0 

2  2  2  2  2  2 

r  3  (x-5)  +  (z-C)  ,  and  R  =  r  +  y  ,  where  J^C^r)  is  the  zero 

order  Bessel  function  of  the  first  kind,  and  the  Integral  sign  ^  is 

0 

to  be  understood  as  the  Integration  along  the  positive  real  axis 
except  for  an  arc  in  the  lower  half  plane  to  avoid  the  positive  real 
root  p  3  a  of  the  denominator* 

From  (4*2)  it  can  be  shown  that 

Gy  -  aG  3  ~  I  in  y  <  0  *  (4*3) 

If  the  integrals  W  and  L  are  defined  by 

W(x,y,z)  =  JJ  f(JLtC)  I  djJdC  ,  (4*4) 

S 

and 

L(x,y,z)  “ JJ  f(Jl»C)H(x,y,z)diidC  1  (4*5) 

S 
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H(x,y,2)  -  J^(pr)«P2rjp  , 

0 

it  follows  from  (4.1)  that  V(x,y,z)  «  W(x,y,a)  ♦  L(x,y,a).  For  a 
function  fCxfS)  continuous  on  S,  the  Integral  L  satisfies  condition 
(3.30),  and  W  +  L  satlsflea  condition  (3»33)»  Furthernora,  by  (4.3) 
one  finds 

V  (x,y,z)  -  aV(x,y,z)  =  W  (x,y,z)  in  y  <  0.  (4.6) 

y  y 

It  is  a  standard  theorem  of  potential  theory  that  W  satisfies  the 
relation 


lim  W  (x,y,z) 

y-*0 

y<0 


f(x,z) 

0 


within  S  y 
outside  S  . 


Hence, 

Vy(x,0,z)  -  aV(x,0,z)  «  f(x,z)  within  S  , 

and 

Vy(x,0,z)  -  nV(x,0,z)  *  0  outside  S  , 


(4.7) 


(4.8) 

(4.9) 


so  that  W  L  also  satisfies  condition  (3*  31)* 

Now  it  can  be  seen  from  (3*32)  that  the  potential  V(x,y,z)  given  by 
(4.1)  will  be  the  solution  of  the  given  boundary  value  problem  if 
f(x,z)  is  chosen  as  a  solution  of  the  following  integral  equation: 


h(xyz)  •  aV(x,0,z)  s  f(x,z) 


(4.10) 


or 
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f(x,z)  JJ  f(]^tC)<>(x,0yZy(,0,C)(i(dC  =  h(x,z)  within  S.  (4*11) 
S 

By  making  use  of  known  integrals  in  [?]«  the  kernel  of  the  Integral 
equation  can  be  explicitly  shown  to  be 

G(x,0,z,(,0,C)  =  ^  -  -naCl^Car)  S^(ar)  -  i2J^(ar)  ],  (4.12) 


where  Y  (ar)  denotes  the  zero  order  Bessel  function  of  the  second 
o 

kind,  and  S^(ar)  denotes  the  zero  order  Struve  function,  respectively. 
Here  one  observes  that  as  ar  becomes  large,  S^(ar)  w  Y^(ar),  hence 
Green's  function  (4.12)  tends  to 


G(x,0,z,j^,0,C)  ~  ■f’ “  2TTa[Y^(ar)  -  iJ^(ar)  ] 

M  —  >  2TTa  ~  ^4)  _  ^  cos(ar  -  tT/4)1 

^  L  /'nar/2  /  TTar/2  J 


(4.13) 


Therefore,  the  kernel  of  the  integral  equation  fluctuates  as  the  fre* 
quency  of  the  oscillation  Increases. 

According  to  the  Fredholm  theory.  Equation  (4.11)  is  soluble 
if  the  homogeneous  equation, 

f®(x,z)  +1^1/  f°(4,C)G(x,0,z,t,0,C)dJidC  »  0  (4.14) 

S 

possesses  only  the  trivial  solution,  f^(x,z)  *  0.  However,  it  has  been 
shown  In  Cl]  that  the  potential  constructed  with  the  solution  of 
Equation  (4.14), 
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V®(x.y.*)  -T^  JJ  f°(4,C)a(x,y.*.5,0,C)dj:dC 


▼anisbcs  identically  in  y  <  0*  Since  V**  vanishes «  one  finds  fron  (4*8) 
that  f^(xys)  must  also  vanish*  Hence,  the  Integral  Equation  (4«ll)  can 

be  solved* 

Upon  finding  the  solution  f(x,z),  then  from  (3>37)  and  (4*10) 
the  pressure  may  be  determined  by 


pg  =  aVg(x,0,z)  a  1  -  fg(x,z),  for  heave 


a  aVj^(x,0,z)  a  z  -  fj^(x,z),  for  roll 

p6  *  aVg(x,0,z)  a  X  -  fg(x,z),  for  pitch 


(4*15) 

(4*16) 

(4*17) 


where  f(x,z)  a  Re[f(x,z)]  *  ImCf(x,z)l*  Then,  substitution  of  (4*15) 
into  (5*54)  yields 


My  a  i  JJ  {l  -  lleCf2(x,z)]}  dS  , 

S 

N  ■  "  ^  JJ  lBi[f2(x,z)]  dS  for  heave. 


(4*18) 


substitution  of  (4*16)  into  (3*33)  yields 

a  -  ^  JJ  z{z  -  Re[fj^(x,z)  ]]■  dS  , 

S 

a  i  JJ  z  Im[fj^(x,z)]  dS  for  roll. 


(4.19) 


i 
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and  subatltutlon  of  (4*17)  Into  (3*36)  ylolds 

*  a  IJ  *{*  “  R«Cfg(x,z)]}  dS  , 

S 

*  “  a  JJ  *  dS  for  pitch, 

S 


(4.20) 


V.  NUMERICAL  PROCEDURE 


Suppose  an  ellipse  * 


2  2 
“If  which  becomes  ■  1 


under  the  trauisformatlon  x  >  i^a,  z  >  z/a,  and  b  *  B/a,  represents  the 
surface  of  a  zero-draft  ship*  A  numerical  procedure  will  be  presented 
which  determines  approximately  the  value  of  the  unknown  density  f(x,z) 
over  the  surface  of  the  ship.  The  procedure  is  based  on  Fredholm's 
method  of  replacing  the  integral  equation  with  a  finite  set  of  linear 
equations  which  relate  the  values  of  f(x,s)  at  chosen  pivotal  points 
on  the  elliptic  surface. 

The  substitution  of  (4.12)  into  (4.11)  yields 
1  b 

2f(x,z)  J  J  f(j:,C)§-tTaCT^(ar)  ♦  S^(ar)  -  i2J^(ar)}dCdt 


-1  -b 


s  2h(x,z)  on  y  ■  S(x,s).  (3*1) 


Since  f(x,z)  is  a  complex  function,  (5.1)  can  be  resolved  into  a  pair 
of  equations 
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1  b 

J  J  log(»r)  H(ar)JdCd^ 

-1  -b 

1  b 

-  J  J  Im[f(  Ji,C)  ]j^(ar)dC<H  ■  2h(x,a)  , 

-1  -b 

1  b 

2ImCf(x,z)  ]  ■•■  a^  J  J  RaCfC^.C)  ]Jp(ar)dCdj; 

-1  -b 


1  b 

*  n  J  J  ^  ti  0  -  a  log(ar)  R(ar)^Cd5  s 

-1  -b 


whera  R(ar)  ■  I  (ar)  +  S  (ar)  -  4  log(ar).  R(ar)  la  a  regular  function 
at  r  >  0* 

A  lattice  la  eatabllahed  on  the  elliptic  aurface  by  dlTldlng 
the  long  axla  Into  alxteen  equal  InterveJLa  h,  and  the  vertical  ordln- 
atea  parallel  to  the  abort  axla  Into  eight  equal  Intervala  k(x)  Caee 
Figure  l]«  Then  each  pivotal  point  can  be  Identified  by  the  coordlnatea 

x^  *  (l-8)h  J  *  0,1,  —  16,  where  h  =  1/8  , 

*j(xi)  ■  (j-4)k(xj^)  j  «  0,1,  —  8,  where  k(Xj^)  »  ^1-x^ 

To  obtain  the  aolutlon  f(Xj^,Zj),  It  la  neceaaary  to  conalder  Equatlona 
(5*2)  only  at  the  forty-one  pivotal  polnta  contained  In  one  quadrant. 
(The  function  becomea  either  aymmetrlc  or  antl-ayametrlc  relative  to  the 
coordinate  axea  depending  upon  the  type  of  oaclllatlon. ) 

Since  r  denotea  the  dlatance  from  a  given  point  P  (x^,Zj), 
where  1  «*  0,1,  —  8,  and  j  «  0,1,  —  4,  to  a  variable  point  Q  (jitO 
*  w4«re  I  «  0,1,  —  16,  and  m  s  0,  1,  —  8,  one  flnda  the 
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Integrands 

f(Ji,C)/y(x^-5)^  ♦  (Zj-C)^,  and  f(ji,C)log 

possess  a  singularity  at  (x. ,z.)  =  (x. ,z  )•  Nevertheless,  these 

1  j  in 

improper  integrals  associated  with  the  above  Integrands  do  exist*  A 
numerical  scheme  will  be  shown  for  evaluating  such  integrals  involving 
two  variables* 


A*  Treatment  of  the  Integral* 


.  f ^  f  -  ■:  awt .  (5.4) 


X  z 
o  o 


y(x^-!i)^+(zj-C)‘ 


In  the  given  lattice,  division  lines  and  curves  form  a  four- 
sided  (two  or  three  sided  near  the  boundary)  mesh*  A  neighboring  region 
6  of  a  point  P  is  defined  as  a  region  which  is  the  union  of  all  meshes 
touching  P*  Then  the  region  S  of  the  Integration  may  be  divided  into 
two  parts,  so  that 


Is(Xi,Zj)  *  I^(x^,z^)  ♦  Ig_^(x^,Zj) 


.5) 


where,  for  example,  if  6  is  a  neighboring  region  about  cm  interior 
point  P, 
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For  th«  ovaluatlon  of  Simpson's  rule  is  gsnsrslly  ussd*  In  csss 

thsrs  are  an  odd  number  of  divisions  in  an  interval  of  integration, 
the  Integral  over  the  division  nearest  to  the  region  6  is  evaluated 
by  the  trapezoidal  rule  instead* 

Now,  setting  t  »  p,  and  C  “  =  q»  (5*5)  nay  be  ex¬ 

pressed  as 


h  k  k 

=  J  J  ^  dqdp  »  J  F(h,q)dq  ,  (5*6) 

-h  -k  ^p^+q^  -k 


IS 

where  g(p,q}  »  f(x.  ♦p,  z.+q),  auid  F(h,q)  «  f  dp*  Since  f(x,z) 

t  J  r?  ? 


2  2 
p  +q 


is  continuous  and  bounded,  it  can  be  shown  that  F(h,q)  is  also  a 
bounded  and  continuous  function  of  q  for  q  /  0*  When  q  approaches  zero, 
F(h,q)  becomes  logarithmically  singular*  Therefore,  a  generalized 
Gaussian  quadrature  formula  is  used  for  an  approximate  evaluation  of  I^* 
To  determine  this  formula,  rewrite  (5*6)  in  the  form 


I 


2 

dq  «  ^  HiF*(h,qi)  * 
i=l 


(5*7) 


The  abscissae  q^  can  be  found  as  the  zeroes  of  the  second  degree  poly- 
nomial  Q(q)  s  1  aq  bq  ,  which  is  a  member  of  the  set  of  polynomials 
orthogonal  in  the  interval  (-k,k)  for  a  weight  function  logiql*  It 
turns  out  that 


<4(q)  «  1 


1  -  31ogk 


q 


2 


f 
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80  that 


q,  =  - 


1  -logk 


and  q. 


,  }sJlzAlp£k 

1-logk 


(5.8) 


Now  the  wolghts  corresponding  to  q^  can  be  determined  by  requiring 
that  formula  (5*7)  be  exact  if  F*(h,q)  is  a  polynomial  of  degree  zero 
or  one.  In  the  present  case  these  are 


=  kClogk-  1)  .  (5.9) 


By  substituting  (5.8)  and  (5.9)  into  (5.7)  one  obtains  a  quadrature 
formula 


Ij(Xj^,Zj)  wk(logk-l) 


J 


g(Pfqi) 


dp  +  J 


^(p.qg) 


L-h  log  q  ^p^+q^  -h  log  ^q^ 


dp 


k(log  k  -  1) 


) 


3C.  , 
4+1 


^  J 


Lxi.i  log  q^{x^~^)^+q^ 


log  q2^*^*i-Ji)^+q|  J 


dC 


(5.10) 


which  gires  exact  results  if  F*(h,q)  is  a  polynomial  of  degree  less  than 
or  equal  to  three* 

Next,  by  Lagrange's  interpolation  formula,  f(x,Zj+q^)  and  f(x,Zj+q2) 
can  be  related  to  the  function  evaluated  at  the  equally  spaced  points 
^(*»*j.l)»  f(x,Zj)  and  f(x,Zj^j^).  Thus,  one  finally  finds  the  for¬ 


mula 
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3^ 

**"  k^log  q-  J  ♦  2(k2-q2)f(Ji,*j) 


"i-l 


♦  q. 


2  2 
+q. 


di: 


(5.11) 


The  integral  expressed  by  (5*11)  is  now  readily  eraluated  by  the 
use  of  Simpson's  rule* 

For  the  boundary  points  (Xj^,s^)  i  *  1»2,  — —  8,  the  associ¬ 
ated  is  treated  in  a  quite  similar  manner*  Here  one  writes 


■,  2 


-  r"  r 


i-l  »o 


dCd^ 


(5*12) 


Following  the  transformations  (5*6)  and  (5*7) i  (5*12)  can  be  expressed 
in  the  form 

h  k  k  2 

VVo’  ■  I  J  ■  I  “  I  V"'*'-''!’' 

-h  oVp%q2  0  i=l 

(5*13) 


Then  the  abscissae  q^  can  be  found  from  the  quadratic  equation 

k  k  k  2  ^  ^ 

[J  log  q  dq  J  q^log  q  dq  -  (J  q  log  q  dq^  Jq^  ♦  J  q  log  q  dq[J  q^log  q  dq 

0  0  0  0  0 

*r  k  k  k  2 

-  J  q^log  q  dq]q  ♦  J  q  log  q  dq  J  q^logqdq-(J  q^log  q  dq)  .0  * 


(5*14) 
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Not*  that  d*p«nds  upon  k(xj^)  In  a  rathar  complicated  manner*  Here* 
the  weights  corresponding  to  are  given  by 


"l  '  <!  'i'!  -  ■’  'io]  •  \ 

I  (5*15) 

®2  *  J^^°®  q  ‘iq  -  log  q  ‘iq]  •  ) 


Hence,  the  quadrature  formula  Involving  the  functions  evaluated  at  the 
equally  spaced  points  f(x,z^}  and  f(x,z^)  becomes 

H-(k-q-) 


♦IrH^(k-q^) 


r\- 


Lk  log  q  j - « — *•  k  log  q.  i - ^ — * 

i-i  ^  ^  V<  V«>  ♦’a 


“i'll 


log  q 


^1-1 


==]f(4,z^)dt.  (5*16) 


The  Integral  associated  with  the  point  (x^z^)  at  the  leading  edge 
of  the  ellipse  also  requires  a  special  treatment*  From  the  definition 
of  the  neighboring  region,  one  starts  by  writing 


2  Q  X..  U  A.  1& 

.'4)  =  r  f  ■ -  J  i  ■ip'iq 


8k  k 


Z  X 

o  o 


2  2 
p  +q 


^  2 

I  -  I  V(p,k)  *  (5*17) 


Then  the  abscissae  p^  caui  be  found  from  (5*14)  by  replacing  the  variable 
q  with  p,  and  the  Integral  limit  k  with  h,  respectively*  It  also  follows 
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from  that  the  weights  corresponding  to  are  given  by 


I<i  ■  p'^'Ip~[p2  J  log  P  dp  -  J  P  log  P  dp]  , 


(5»18) 


^2  “  P5-P 


aa  aa 

““[pj^  J  log  p  dp  -  J  p  log  p  dp]  . 


The  quadrature  formula  can  now  be  written  as 


*p8^Pl^  I^(h-p^)  f(x^,C)  *p8^^2^  L2(h-P2)  f(x^.C) 

J  h  log  p^  ""ry  -  J  h  log  p-  Ip' 


log  Pp  1-5 - -J 


r®  LiPi  f(xj^,c)  ^2  LgPa  f(xj^«C) 

■*■  J  hlogp^  - 2  •'  hlogP2  fY~ - 72 


(5.19) 


Since  Z|^  a  0,  and  since  f  [x^,Zj(pj^)  ]  =  f(x^,Zj^)  j  ■  0»1»  8t 
(5*19)  may  be  rewritten  as 


r*p8^Pl^  ^(h-p, )  ,  *p8'P2^  Mh-pp)  ,  1 

V*o»%^'^ir  h  log  p  h  log  P2  -r:==dCjf(x^.z^) 

‘-Zo(Pi)  *o<P2>  ^ 


*p8^P2^  L-(h-p 


J  hiogp  rr— 2  ^  *  J  hiogpp  /  2  2 


(5.20) 


Here,  functions  f  [xj^,Zj(pj^)  ]  may  be  approximated  by  f  [xj^,Zj(Xj^)  ]  ^or 
related  to  fCx^pZ^Cx^)]  by  linear  interpolation^  to  evaluate  the  third 


and  fourth  integrals* 
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B*  Treatment  of  the  Integral* 


\  f(4,C)iogy(x^-t)  +(zj-C)^  dCdj; 


X  z 
0  o 


(5*21) 


If  the  region  S  of  integration  Is  divided  into  a  neighboring 
region  6  and  a  remaining  region,  one  obtains 


Jg(Xi,Zj)  =  +  ’^s-6^*i»®j^ 


Applying  Simpson's  and  the  trapezoidal  rules,  the  Integral  *^3.^  can 
be  evaluated  by  a  process  similar  to  the  one  used  for  the  evaluation 
of  ^3.5*  For  example.  If  6  Is  a  neighboring  region  about  interior 
points  P,  the  integral  is  given  by 

^+1  *J+1 

J^(Xj^,Zj)  =  J  J  f(ji,C)logy(x^-ii)^+(Zj-C)^  dCd^  .  (5.22) 

*1-1  *d-l 


By  the  transformation  j^-x^  =  p,  C-*j  ■  q,  (5*22)  may  be  rewritten  as 

h  k  _  k 

*  J  J  s(p*<l)logyp^+q^  dpdq  =  J  a(h,q)dq  ,  (5»23) 


-h  -k 


-k 


where  g(p,q)  «  f(x^*p,  ^j+qJ*  Q(h,q)  =  J  g(Ptq)log^^H^  dp  .  It 
can  be  shown  that  Q(h,q}  Is  not  only  a  bounded  and  continuous  function 
of  q  for  q  ^  0,  but  also  remains  finite  when  q  approaches  zero*  There¬ 
fore,  non-generallzed  Qausslan  quadrature  formula  may  be  used  to  evalu- 
ate  Jj*  In  order  to  derive  such  a  formula,  write  (5*23)  in  the  form 


(5*24) 
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fV 

■  J  a(h,q)dq  m  ^  M^G(h,q^)  . 

-k  i=l 


Th«  abscissa*  q^  are  obtained  as  seross  of  th*  second  degree  poly¬ 
nomial  Q(q)  s  1  *  q^;  thus  q,  »  -  — ,  and  q,  *  —  .  (5.25) 

^  ^  /5 

Now,  the  weights  corresponding  to  q^  are 


Ml  »  =  H 


(5.26) 


Then,  substitution  of  (5.25)  and  (5.26)  yields  a  quadrature  formula 
'^6^*1*® lt[J^K(p»qi)loS^^*^^  dp  ♦  j^g(p.q2)lo«E^^^  dp] 

»  k[j  f(5,Zj+qj^)logy(Xj^-5)^+q2  dj^  +J 


‘i-1 


‘i-1 


(5.27) 


which  again  gives  exact  results  if  G(h,q)  is  a  polynomial  of  degree 
less  than  or  equal  to  three. 

If  f(x,Zj-*'q^)  and  are  expressed  in  terms  of  f(x,Zj^^), 

f(x,Zj)  and  f(x,Zj^^)  by  Lagrange's  interpolation  formula,  one  finally 
finds 

^+1 

«  J  J  [q2^(ll.*j.i)+2(k^-q2)f(i;,Zj)+q2f(i;,Zj^j^)]l08y(x^-t)^+q|  d^ 

^i-1 

'  f  1  [fUi*j.i)+^f(4i®j)-''f(Ji**j+l)]lo8y(x^|-I)^^  d^  .  (5.28) 

^i-1 
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Not*  that  quadrature  coefficients  are  equal  to  those  of  Simpson's  rule* 
For  boundary  points  (xj,s^)  i  =  i»2»  "•*“  8»  the  integral  ia 

J^CXi^z*)  »  f  f(i,OlogJ(x^^0^*(\~C^  dCdil  •  (5.29) 

*i-l  *0 


Then  (5.29)  can  be  transformed  into 


r—  ^ 

J.(x.  ,z  )  =  I*  /'g(p,q)fnX^r?  dqdp  »  rG(h,q)  «  ^M  G(h,q^).  (5*30) 

*  ^  ®  ih  •'o  0 


In  the  present  cassi  the  abscissae  q^  and  corresponding  are  readily 
found  as 

and 

=  H2  =  k/2  . 


(5.31) 


Hence,  in  terms  of  f(x,z^)  and  f(x,z^), 
can  be  escplicitly  shown  as 

*i-l 


the  quadrature  formula  (5.30) 


^]f(ji,z^)dli 


X  _ 

+  losy(x^-Ji)^+q^  +-^  *  (5.32) 

*i-l 


Note  that  quadrature  coefficients  are  similar  to  those  of  the  trape** 
zoidal  rule.  The  Integral  associated  with  the  boundary  point  (x^,Zj^) 
requires  a  special  treatment*  Now,  a  quadrature  formula  will  be  derived 


by  writing 
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0  0 


*  J  dj^dC  =  f  J* g(p.q)log7p^+q^  dpdq 

^  ^  z  X  ®  ^  0  0 

0  o 

U  2 

=  j  G(p,k)dp  «  ^N^a(pj^,k)  .  (5.33) 


(5.33) 


Th«n,  from  (5.31)  one  finds 


(5.34) 


-  h/2  . 


Since  f[x^,*j(p^)]  s  f(x^,zj^)  j  s  0,1,  —  8,  one  obtains  the  quadra¬ 


ture  formula 


'*4’“  f  — ji-logJpJ-C^  dC*J*  dC  t<»  ,1 


r^^l^  p  p— —  p  r-“~” 

“T^os/pf+C^  f(x.,C)dC+J  -SlogJ^^  f(x.  ,C)dC  . 
<Pi)  ^  ^  ^  ^  igCp,)  ^  ^ 


(5.35) 


S***®*  functions  f  [xj^,Zj(Pj^)  ]  nay  be  approximated  by  f  [xj^,2j(xj^)  ]  (or 
related  to  f[xj^,aj(xj^)  ]  by  linear  interpolation.) 

Application  of  the  singular  integral  formulae: 


for  inner  points 
for  boundary  points 


(5.11)  and  (5.28), 
(5. 16)  and  (5.32), 


for  leading-edge  points  (5.20)  and  (5.35)» 
together  with  Simpson's  and  the  trapezoidal  rules  enable  one  to  reduce 
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th*  pair  of  integral  aquations  (5»2)  to  alghty-two  linear  equations 
relating  the  values  of  SeCfCxts)]  and  ImCfCxyz)]  at  forty-one  pivotal 
points* 

For  example,  if  (x^,*j)  are  inner  points,  (5*2)  is  reduced  to  the 


form 


i+1  j+1  '' 

l«i-l  m=j-l  maO  J=i+1  maO 


i-1  8  16  8 


i+1  J-1  i+1  8  I 

•I  I*  I  I  "iJ 


l«i-l  maO  Ni-l  m=j+r 

)=0  maO  2aO  maO 

16  8  I  i+1  j+1 

J»0  m=0  iai-1  m=J-l 

•  i-1  8  16  8  i+1  j-1  i+l  8 


im 


pi-i  o  lb  o  1+1  J-1  1+1  o  T 

[il*I  1*1  i*  1  I  los  "l  J 

LlaO  maO  2=i+l  maO  iai-1  maO  }ai-l  maj+l-* 


16  8 

'I  I^2m  ^  R(arj^)[lmCf^“]  .  0 

2a0  ni=0 


(5*36) 


where  the  notations  represent,  respectively, 

f^^  a  f(x^.Zj).  f^®  a  f(Xj,Z^),  r^^  a  X^-Xj )^  +  ( Z j-Z^)^ 

*‘2m  •  ^/^  *T  TT~Hof  kj  ’  '  h(x^,Zj),( 
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h.  ■  h  D^k(x^)  (Th*  coefficients  for  Simpson's  or  the 

trapesoldal  rule). 


(5.37) 


2  2  Z 

c  h  ^  or  C  s  "  .  C.h  D  k(x  ). 

I  k  log  qg  i  k  log  q2  *  )  m  I  m  I 


In  (5*36),  am  appropriate  change  Is  necessary  If  (x^,Sj)  are  boundary 
points*  Furthermore,  one  observes  that  the  function  f(x,s)  has  the 
following  properties; 

for  h.*<.  _  \ 

for  roll  f*."  .  I  (5.58) 

for  pitch  f*’“  •  j 

Therefore,  using  (5*33)  the  left  hamd  side  of  (5*36)  can  be  formed  Into 
a  rectangular  matrix  for  unknowns  f^^*  In  this  mamner,  one  obtains 
three  sets  of  linear  equations  describing  the  type  of  forced  oscillation. 
For  vamious  values  of  the  frequency  parameter  a  «  ao  /g  ■  2Tiai/\,  one 
cam  determine  the  values  of  the  density  function  by  solving  these  equa¬ 
tions* 

Actual  steps  of  the  computation  work  for  obtailning  the 
solution  consist  of: 

[step  l]  Determination  of  the  distamce  r  from  the  given  pivotal  point 
to  all  pivotal  points  on  the  lattice* 

[step  2]  Cailculatlon  of  functions  appeairlng  In  the  coefficients  of 

the  linear  equations  by  seams  of  Interpolation  from  the  given 
table* 
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[St«p  3]  Sumaatlon  of  the  coefficients  and  grouping  of  the  matrix 

In  accordance  with  the  type  of  oscillation* 

[step  4]  Numerical  solution  of  linear  equations  by  the  elimination 

process  based  on  the  algorithm  of  Gauss* 

VI.  NUMERICAL  RESULTS 

The  results  presented  in  this  paper  are  obtained  by  the 
application  of  a  new  numerical  scheme  developed  for  evaluating 
improper  double  Integrals  associated  with  the  fundamental  singular¬ 
ities  1/r  and  log  r .  As  it  will  be  shown,  the  present  results  agree 
I*lrly  well  with  the  results  computed  by  the  application  of  the  pre¬ 
vious  scheme  (see  Section  V  in  [8])  which  uses  a  circular  neighboring 
region  about  the  singularity*  Secondly,  in  the  present  work,  the 
added  moment  of  inertia  and  damping  coefficient  for  roll  and  pitch  are 
computed  correctly*  These  quantities  presented  in  [8]  are  in  error* 

It  was  discovered  that  a  mistake  was  committed  in  Step  3  of  solving 
linear  equations  for  the  source  density*  In  the  course  of  grouping 
the  matrix  in  accordance  with  rolling  or  pitching  oscillation,  the 
signs  securing  in  (3*38)  were  unwittingly  reversed*  Subsequently,  the 
added  moment  of  Inertia  euid  damping  coefficient  were  determined  from 
erroneous  values  of  the  pressure*  Lastly,  in  the  present  work,  the 
range  of  the  computation  was  extended  up  to  the  frequency  parameter 
a  «  2rr.  However,  as  has  been  seen  in  (4*13)  because  of  the  fluctuating 
tendency  of  the  kernel  at  large  values  of  the  argument,  the  accuracy 
of  results  decreases  beyond  a  >  tr.  To  assure  the  same  accuracy  over 
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th«  range  of  a  >  to  2tt,  the  mesh  spacing  may  be  bisected* 

The  investigation  was  performed  for  am  elongated  ellipse  of 
axes  ratio  b  s  b/a  «  1/8 «  and  1/4,  and  its  limiting  case,  a  circle  of 
b  B  1.  Thus,  three  groups  of  computations  were  necessary  with  each 
value  of  the  frequency  parameter  a,  which  assume  the  values  w/6,  tt/^, 
V4,  Ti/3,  2tt/5»  tt/2,  2n/3,  tr,  5tt/4,  3tt/2,  7V4,  and  2tt. 

For  each  combination  of  b  and  a,  three  sets  of  linear  equa¬ 
tions  describing  the  type  of  oscillations;  namely,  heave,  roll  and 
pitch,  were  solved*  From  the  solutions  of  the  linear  equations, 
pivotal  values  of  the  pressure  were  determined*  Then,  using  Simpson's 
rule  one  finds  normalized  added  mass  M,  normalized  added  moment  of 
inertia  I,  and  normalized  damping  coefficients  N  or  H  by  the  pair  of 
appropriate  formulas  (4*18)  to  (4*20)* 

In  order  to  substantiate  the  preceding  remark,  the  values 
of  M  and  N  computed  by  the  use  of  the  previous  scheme  and  those  com- 
puted  by  use  of  the  present  scheme  are  compared  in  Table  1* 

TABLE  1 


Previous  Result 
b«l  bal/8 


a 

My 

Ny 

My 

Ny 

V5 

2*129 

1*016 

0.102 

0*032 

tt/4 

1*987 

1.026 

0*098 

0*036 

V3 

1*809 

1*002 

0*090 

0*041 

2tt/5 

1.706 

0*968 

0*085 

0*043 

V2 

1.595 

0*910 

0*077 

0*045 

2tt/3 

1.473 

0*811 

0*066 

0*044 

TT 

1.341 

0*628 

0*051 

0.039 

ba 

Present 

1 

Result 

b= 

1/8 

My 

Ny 

My 

Ny 

2*118 

1.022 

0*113 

0*031 

1*981 

1.031 

0.108 

0*034 

1.807 

1.006 

0*099 

0*038 

1.704 

0*970 

0*093 

0*040 

1.592 

0*909 

0*084 

0*041 

1.467 

0*807 

0*072 

0*040 

1.319 

0*658 

0.056 

0*033 
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In  Figur#  2  and  Figure  3«  the  dependence  of  — j—  M  and 

-2  ^  y 

N  on  the  parameter  a  are  presented*  In  Figure  4  and  Figure  5f 
the  quantities  I  and  H  are  plotted  as  functions  of  the  para- 

^  X  AX  ^2  ^2 

TTA 

met«r  a*  Similarly,  the  quantities  plotted  In 

Figure  6  and  Figure  7%  respectively.  The  multiplication  factor 

a 

“X"  ®  ordinates  represents  the  ratio  of  the  area  of  a  circle 

^  D 

with  radius  equal  to  a  (half  the  long  axis  of  the  disk)  to  the  area  of 
the  waterplane  of  the  disk  under  consideration.  This  factor  was  intro¬ 


duced  In  order  to  make  the  curves  comparable.  The  actual  vad-ues  corres¬ 
pond  to  one-fourth  of  the  ordinate  for  the  ellipse  of  b  =  1/4,  and  one- 
eighth  of  the  ordinate  for  the  ellipse  of  b  =  1/8,  respectively. 

The  curves  for  the  circular  disk  b  =  1  in  Figure  2  and 
Figure  3  compare  very  closely  to  the  corresponding  curves  in  Figure  6 
and  Figure  7  of  C5],  which  were  obtained  by  treating  the  circular  disk 
problem  by  deleting  the  angular  dependence.  The  same  method  Is  applied 
to  solve  the  case  of  pitching  in  [61.  In  Table  2,  the  results  obtained 
by  the  present  method  and  those  obtained  by  the  method  of  [61  are  com¬ 
pared.  In  the  latter  method,  the  accuracy  of  the  computation  Is 
Increased  by  meeuxs  of  doubling  the  numbers  of  pivotal  points  on  the 
radius  of  the  circular  disk.  Over  the  range  of  a  =  tt  to  2tt,  the 

present  method  gives  the  lower  values  of  I  and  H  which  amounts  to 

z  z 

3  to  14^  for  I^  and  8  to  20^6  for  H^,  respectively. 
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TABLE  2 


For  a  Pitching  Circular  Disk 


By 

Present 

Method 

By  Method  of  C6] 

a 

I 

H 

I 

H 

z 

z 

z 

z 

tt/6 

0*281 

0*013 

0*282 

0*014 

n/5 

0*285 

0*020 

0*287 

0*021 

tt/4 

0*289 

0*031 

0*293 

0*034 

tt/3 

0*289 

0*051 

0*294 

0*055 

2n/5 

0*284 

0*064 

0*290 

0*070 

it/2 

0*273 

0*079 

0*279 

0*085 

2V3 

0*252 

0*091 

0*259 

0*097 

TT 

0*219 

0*091 

0*227 

0*099 

5V4 

0*203 

0*086 

0*212 

0*095 

3tt/2 

0*190 

0*079 

0*203 

0*091 

7Tt/4 

0*176 

0*071 

0*198 

0*085 

2tt 

0*167 

0*063 

0*194 

0*079 

It  should  bs  noted  that  as  the  frequency  of  the  forced  oscillation  o 

2 

tends  to  aero  (hence,  a  **  0),  M  becomes  a  constant  and  N  (being  a  S  ) 
the  damping  factor  will  vanish*  Since  M  is  0(log  a),  in  the  two- 

y 

dimensional  case,  M  tends  to  infinity  while  N  becomes  a  constant  as 

y  y 

the  frequency  tends  to  zero.  In  Figure  3,  the  limiting  values  of  the 
ordinate  may  be  shown  by  the  relation 


lim 

a-*0  A 


N 


where  A  ■  rmS,  and  S 


Tia^  /  -2.  b  2 

(a  S  }  ■  —  TT  a 

a 


^  2 

TwiB/a*^  ■  TiB/a  for  the  elliptic  disk* 


Hence, 


at  the  origin,  the  slopes  are  given  by 
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for  b  »  1  TT^  -  9»8696  , 

for  b  -  lA  IAtt^  =  2.4674  , 

for  b  »  1/8  1/8 TT^  »  1.2357,  roopoctively. 

Th«  computed  results  dsTlate  very  rapidly  from  the  low  frequency  approzi* 
matlon  shown  in  the  above.  Furthermore,  as  the  frequency  tends  to  zero 
in  Figure  4  and  Figure  6,  amd  I  become  constants,  while  in  Figure  5 

X  z 

and  Figure  7,  H  and  H  will  vanish. 

X  z 


This  work  was  performed  at  the  Boeing  Scientific  Research  Laboratories 
based  on  the  research  initiated  with  the  support  of  the  Office  of  Naval 
Research  Contract  NONR  760(21)  at  Carnegie  Institute  of  Technology. 

The  author  expresses  sincere  appreciation  to  Professor  R.  C.  HacCamy 
for  his  valuable  advice  and  criticisms,  and  to  Professor  T.  E.  Stelson 
for  his  support  and  suggestions.  The  author  is  also  Indebted  to 
Dr.  J.  F.  Price,  who  assisted  for  developing  new  numerical  procedure; 
to  Dr.  T.  E.  Turner,  who  has  supervised  the  continuing  research  at 
the  Boeing  Scientific  Research  Laboratories;  and  to  L.  N.  Hakala,  who 
programmed  the  computation  work  for  the  IBM  7090  data  processing  system 
at  the  Boeing  Company  in  Seattle,  Washington. 


39 


BIBLIOGRAPHY 


1*  Stoker  and  Peters,  "The  Motion  of  a  Ship  as  a  Floating  Rigid  Body 

In  a  Seaway,"  Conunun*  on  Pure  and  Applied  Math*,  Vol*  X,  No*  3,  1937* 

2*  John,  F* ,  "On  the  Motion  of  Floating  Bodies,"  II*  Coaunun*  on  Pure 
and  Applied  Math.,  Vol*  III,  No*  1,  March  1950* 

3*  Hasklnd,  M*C*,  "The  Rydrodynamlcal  Theory  of  the  Oscillation  of  a 
Ship  In  Waves,"  Russian  translation,  1946* 

4*  MacCamy,  R*C*,  "On  the  Scattering  of  Water  Waves  by  a  Circular 

Disk,"  Archive  for  Rational  Mechanics  and  Analysis,  Vol*  3,  No*  2, 

1961* 

5*  MacCamy,  R*C*,  "On  the  Heaving  Motion  of  Cylinders  of  Shallow  Draft," 
Journal  of  Ship  Research,  Vol*  5,  No*  3,  1961* 

6*  Kim,  W*D*,  "The  Pitching  Motion  of  the  Circular  Disk,"  (to  be  published)* 

7*  Grobner  and  Hofrelter,  "Intergraltafeln,"  Sprlnger-Verlag,  3rd  Ed*, 

1961* 

8*  Kim,  W*D*,  MacCamy,  R*C*,  and  Stelson,  T*E*,  "The  Forced  Oscillation 
of  Shallow  Draft  Ships*"  Final  Rep*,  Contract  No*  N0NR-760(2l) , 

Carnegie  Institute  of  Technology,  Pennsylvania  (1962)* 


nOURE  3 


5 


1 


3 


6 


7 


